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ABSTRACT 


Diagnostic model output parameters, provided by the Fleet 
Numerical Weather Central, Monterey, California (FNWC), and 
the marine fog frequency climatology developed at the Naval 
Postgraduate School, Monterey, California, are statistically 
processed in context with marine surface synoptic ship re- 
ports for the purpose of developing a linear regression 
scheme for modeling the distribution of marine fog. The study 
area encompasses a large section of the North Pacific Ocean 
(from 30-60N). The time period involves 0000 GMT, 1-30 July 
1976. The predictand is a fog parameter developed from a 
quantitative categorization of each of the 4481 synoptic ob- 
servations according to reported present and past weather, 
low-cloud type, and visibility codes. The 38 model output 
and climatological predictor parameters are interpolated to 
each of the ship observation positions and the resultant data 
file is used for the derivation of the regression equations. 
The diagnostic capabilities of the regression equations, 
along with other existent approaches, are analyzed through the 
use of three verification scoring systems--Heidke Skill, 
Threat, and the Panofsky-Brier Probability scores. Improve- 
ment over climatology and FNWC's operational fog probability 
program (FTER, applied at analysis time) is demonstrated. 
selective mappings of the regression equation outputs and 
categorized observations are intercompared with the sea-level 
pressure analysis; FTER; and the most significant predictor 


parameter in the regression equations, evaporative heat flux. 
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I. INTRODUCTION AND BACKGROUND 


Marine fog continues to be a threat to safe nautical and 
lov-level aviation activities, civilian and military. Vheeler 
and Leipper (1974) documented various monetary and human 
losses associated vith United States Navy operations in poor 
visibilities due solely to fog. An accurate analysis of 
marine fog, and its associated increases in forecasting 
accuracy, should help decrease these losses as vell as en- 
hance the specifications of lov-level cloud inputs to the 
various planetary boundary layer models in the environmental 
sciences community. 

Research on marine fog by the Departments of Meteorology 
and Oceanography at the Naval Postgraduate School (NPS) has 
concentrated on the climatology and the diagnosis of this 
phenomenon from observational data (Renard, Englebretson, and 
Maughenbaugh, 1975; Renard, 1976). A recent effort (Ihli and 
Renard, 1977) attempted to utilize satellite data in the 
diagnosis of marine fog regimes. The results indicate that 
the approach has potential, but more data are needed for fur- 
ther testing. In the realm of forecasting marine fog, an 
earlier study (Nelson, 1972) attempted an approach similar to 
the Model Output Statistics (MOS) method developed by the 
National Weather Service (NWS) (Glahn and Lowry, 1972). 
Nelson correlated visibility at sea over the North Atlantic 


Ocean to locally observed parameter values and developed 


11 





regression equations using the "perfect prognostic" approach. 
Out of the 20 or so parameters tested, vind speed, relative 
moisture content,and parameters associated vith evapora- 

tion shoved the highest correlations vith visibility. He then 
tried to forecast visibility using, as predictors, grid point 
output values of these parameters from the operational models 
at the United States Navy Fleet Numerical Weather Central, 
Monterey, California (FNWC). The results showed little skill 
in the forecast mode. 

Only FNWC is presently producing fog forecasts on a hemi- 
amerite basis through their statistical probabilistic product 
called FTER (see Appendix A). This product is the result of 
calculations based on a statistical processing of certain 
model output parameters related to fog occurrence (U. S. Naval 
Weather Service Command, 1975). FTER does not use any type 
of climatological parameter, partly due to the unavailability 
Of accurate climatologies at the time of the development of 
mae product. Only a limited evaluation of the FTER product 


has been accomplished to date (Renard, 1975b). 
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EE e RN APPROACH 


The primary obiective of this study is to utilize the MOS 
concept in the development of a statistical equation (or set 
Of equations) that vould accurately specify the distribution 
o£ marine fog over an extensive oceanic area, such as the 
North Pacific Ocean. The secondary objective is to evaluate 
the skill of an existing operational scheme, FNWC's FTER 
product, and the NPS £og frequency climatology for the same 
region (Englebretson, 1974, Daughenbaugh, 1975, Villms, 1975). 
In essence, this is an "imperfect diagnostic" approach de- 
veloped in part from numerically analyzed/predicted atmos- 
pheric and oceanographic parameters, similar to the MOS 
concept. 

The approach used in pursuit of the primary objective 
involved three steps: a) the classification of each synoptic 
observation in the data base (described in Section IV) by 
cloud, weather, and visibility reported (predictand); b) the 
interpolation of model output grid-point parameters (and 
other parameters derived therefrom) and a climatological 
parameter (predictors) to the observation points; and c) the 
search for statistical relations between these predictors 
and the predictand. The immediate goal was to find one or 
more multiple, linear regression equations that would be 
Significantly more skillful than either the NPS climatology 
or the FNVC FTER product in describing the actual distribu- 


tion or fog at analysis time. 
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HİT. DATA 


Er AREA OF STUDY 

The NPS climatological fog frequency study (Willms, 1975) 
was confined to the North Pacific Ocean. Since this was 
judged to be the most accurate climatology of fog occurrence 
to date, the same general area was chosen for this study. 
Further restrictions were dictated by the desire to avoid 
the influences on fog patterns by large land masses and the 
large-scale upwelling in the waters off of the western coast- 
line of the United States. The area of study is displayed in 
Fig. l on a Mercator projection. The area is outlined using 
the boundary I,J points of the standard 63-by-63 grid used 
by FNWC. Note the resultant curvature of the grid lines when 
the grid is placed on a Mercator projection. Fig. 2 has the 
same area and outline displayed on a Polar Stereographic pro- 
jection. The entire FNWC I,J grid, on a Polar Stereograpnic 


projection, is shown in Fig. 3 for reference. 


17 Lie PERIOD FOR STUDY 

The month of July, 1976 was chosen as the basic time 
period. The main reason for this selection is the high 
relative frequency of marine fog for this month, as displayed 
in the NPS climatology (Willms, 1975). 

The 0000 GMT synoptic time was chosen as best since this 
time is local noon over the international dateline (180 de- 


grees) vith daylight hours occurring over all of the area of 
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study. This factor should give maximum relative accuracy of 
observations by the transient ship observers (Nelson, 1972). 

Other aspects involved in using only 0000 GMT data are 
persistence and diurnal variations. It is widely held that 
diurnal variations of meteorological phenomena in oceanic 
regions are minimal. In any case, the use of one synoptic 
time will eliminate any diurnal effect. 

Panofsky and Brier (1968) suggest that, in order to 
eliminate persistence at stationary observation points, data 
from every third or fourth day be used in the dependent 
sample. Since the vast majority of the observations used for 
this study were from transient (non-stationary) ships, such 


persistence effects are considered negligible. 


eee oY NOPITIC WEATHER REPORTS 

FNWC provided all of the primary-time synoptic weather 
observations from their data archives. Duplicate reports and 
reports obviously in error were eliminated from the observa- 
tions provided. In all, 4481 0000 GMT observations were ob- 
mained from the period 1-30 July 1976. Observations for 
31 July 1976 were not available from FNWC. 

While the majority of the observations came from tran- 
Sient ships, Ocean-Weather-Station P and ll land (island or 
coastal) stations (see Fig. 1) were also included in the 4481 
observations. The ll land stations are listed in Appendix B, 
Ewen with its station number, elevation, and location (lati- 


26 longitude). The location of land station number 32618 


15 


ә 
a 





was adjusted from its actual location (see Fig. 1) to just 


inside the study area. 


57 MODEL OUTPUT PREDICTOR PARAMETERS 

FHWC also provided all of the selected model output para- 
meters (MOP's) for the time periods involved. These MOP's, 
their description and/or applicable remarks, and range of 
values (from histograms not shown) are listed in Appendix C. 
Some o£ the MOP's were chosen after consultation with Mr. Leo 
Clarke, Director of Research, FNWC, and others were chosen 
using the experience and advice of earlier studies (Schramm, 
1966; Welson, 1972; Grisham, 1973, Renard 197/5a). The only 
desired MOP not available from the FNWC archives is a dew- 
point variable, but the EAIR (vapor pressure) parameter used 
in the study is derived from the dew point. The acronyms 
or abbreviations listed in Appendix C will be used as the 
variable names or symbols throughout the ub. of the 


report. 


Pee DERIVED PREDICTOR PARAMETERS 
Other parameters used in the study were derived from the 


- These are listed 


basic set of 21 MOP's provided by FNWC. 
in Appendix D in the same manner as the basic MOP set is 
listed in Appendix C. As with the basic set, the acronyms 


25 P-bbreviations listed are used in the remainder of the 


report. 


cent Sensible Heat Flux (SHF) which was recovered 
from SHF = SEHF - EHF; it is listed in Appendix C. 
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F. CLIMATOLOGICAL FOG FREQUENCY PREDICTOR 

The values of the NPS July climatological fog frequencies 
(Willms, 1975) used in the study were determined from inter- 
polations to a 249-by-249 grid, which is about 1/4 of the 
FNWC I,J grid (Fig. 4). The resolution afforded by this 
fine spacing enabled the researcher to define all of the major 
Features on the NPS July climatology charts. Such spacing is 
compatible with the original development of the fog fre- 
quencies, which were organized by one-degree-by-one-degree 
grid squares. Some subjective extensions to Willms' (1975) 
final isopleths were necessary in order to accommodate to 
the study area. These extensions are shown as the dot-dash 
meopleths in Fig. 4. 

The overlay used for determining the climatological values 
at the grid points of the 249-by-249 grid was developed by 
computer plotting latitude, longitude values from a modified 
FNWC 1,J-to-latitude/longitude program (See Appendix Ei, 

The other predictor variables were available for the standard 
63-by-63 FNWC I,J grid-point locations. The abbreviation, 
CLIMO, used in the following portions of the text, stands for 


the NPS fog frequency climatology. 
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IV. POGCAT PROGRAM 


A. BACEGROUND 

Until a methodology for using satellite data to identify 
marine fog regimes becomes operational, the main information 
sources for this purpose are synoptic ship reports and some 
island station reports. Previous researchers at the NPS used 
three elements of a ship"s synoptic veather observations to 
categorize or differentiate betveen reports indicating fog 
("foggers") and those indicating veather other than fog ("non- 
foggers"). The three elements used vere present veather (vv), 
past weather (W), and visibility (VIS), See Tables I and II. 
While these elements certainly can be used to ascertain the 
occurrence/non-occurrence of fog, at least one additional 
element appears useful in determining the existence of fog 
regimes, namely, low cloud type (CL). See Table I. 

Some of the combinations of ww, W, and CL used here as 
indicators of the occurrence of marine fog generally would 
not be considered valid in most verification schemes. Hov- 
ever, this non-verification is primarily due to the existing 
rules in veather observing (U. S. Department of Commerce, 
1969). These rules require the reporting of precipitation 
(ww more than 49, such as light drizzle or rain), even though 
a co-existing non-precipitation phenomenon (such as one of 


the modes of fog in the ww 40s group) may have more operational 


Significance. 
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B. CATEGORIZATION PROGRAM 

A scheme (FOGCAT; Fog Categorization) was devised to place 
each of the 4481 observations into one of the five major cate- 
gories: Strong foggers (S), Foggers (F), Past/veak foggers 
(P), Maybe foggers 010), or Non-foggers (N). The scheme 
evaluates only the present weather (ww), past weather (W), 
and low cloud type (CL) codes for category placement. The 
synoptic visibility (VIS) code is not considered until the 
score assignment phase. A description of the FOGCAT scheme 
follovs. 

The 72 combinations of vv, W, and CL that are used in 
categorizing the observations are listed in Table IV, using 
the legend listed in Table III. The partitioning of the 
major categories into subcategories vas done for three pur- 
poses, the least important of vhich is the collection of 
occurrence statistics for each combination. The listing order 
of the subcategories in Table IV is a subjective estimate o£ 
the degree of certainty of a given combination indicating the 
presence of operationally important fog at an observation 
point. The listing order was arranged with a view to both 
u ... or thc observer and the precedence rules for report- 
ins. In effect, it is assumed that the more and/or the 
denser the fog indicators reported by an observer (trained or 
2 .1ned) in an observation, the more certain it is that the 
observation was taken in or near an existent fog regime. The 
second aspect of the listing order is that it provides a 


means of assigning a probabilistic score to each such 
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epeervarion Of ww, W, and CL; that is, a numerical value of 
the degree of certainty that fog was occurring. 

The main problem encountered in devising FOGCAT involved 
the grouping of the various values of ww, W, and CL. The 
groupings used are listed in Table III, which is also the 
legend for Table IV. There are several major differences 
between the present (ww) and past (W) weather groupings 
therein and those used in the development of the NPS fog fre- 
@@enecy climatology (Willms, 1975). The 50's (drizzle) and 
some of the 60's (rain) code values (see Table I) for ww are 
separated from the "any others" group. The ww code values 
of 10 (deep light fog), 26 (heavy fog in past hour), and 40 
(light fog) are placed in the same group (106). The code 
values of 20 (drizzle in past hour) and 24 (freezing drizzle 
in past hour) for vv are moved from the "any others" group to 
the one (11G) containing the values of 11 and 12 (shallov 
heavy fog). The elements in each of the latter tvo groupings 
do not necessarily have similarities in time or nature of fog 
but are related herein primarily due to a subiective consid- 
eration of their respective indicative qualities for the 
diagnosing of fog regimes. 

The W code value of 5 (drizzle) is paired with 4 (fog) 
and this group (4,5 in Table 111) is treated as a reliable 
fog indicator only vhen in combination vith the vv groups of 
41G (deep heavy fog) or 10G (the latter only for major cate- 
gory S). Othervise, the V value of 5 is grouped (5G) with W 


values of 0-2 (sky coverage changes), 
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The lov cloud (CL) type groupings are more straightfor- 
ward. Values 5 (stratocumulus) and 7 (fractostratus of bad 
weather) are used to counter the possibility of an untrained 
observer substituting either of these for the more properly 
reported 6 (stratus or fractostratus). CL values reported 
mem, ) indicate am obscuration. Of the 4481 observations 
used in the study, 919 have a CL of (/), and 466 of these are 
in the a2 subcategory. 

1. Major Category 5 

The ww, W, and CL combinations comprising the strong 
fogger category generally have at least two strong fog indi- 
cators in combination. The strongest combination indicating 
fog is given by subcategory al. However, over 50 per cent of 
the total S cases, 793, are categorized in the aZ subcate- 
gory. A relatively high number of all foggers vere expected 
to fall into this subcategory because of the obscuring nature 
of fog and the results of Willms' (1975) study. 

2. Major Category F 

in general, there are at most two fog indicators in 
the combinations comprising this category. Placing subcate- 
gory dl into the F, rather than S, category indicates the 
20 0£ s opinion that the CL value of Bis not as certain a 
fog indicator as the CL value of 6 in subcategory cl, Sub- 
categories d3 and d4, which are the first occurrences in the 
listing order of a 50G ww value, can be considered equivalent 
to al and a2 in a fog situation except that the precedence 


reporting rule prohibits reporting one of the 41G values. 
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But the total number of such cases is small and the effect 
on the final regression equation should also be minor. 
3. Major Category P 
The past/weak fogger category consists mainly of com- 
binations which have a W value of 4 and no strong indication 
of fog at or within the past hour of observation tine. 
4. Major Category M 
The "maybe" category represents the gray area of the 
scheme. Some of the reports that are placed in this category 
are likely to be in or near a fog regime, but are not veri- 
fiable as fog. 
5. Major Category N 
The non-fogger category was expected to include the 
majority of all the cases, and it did, as evidenced by 2004 


observations or 56 per cent of the total. 


C. SCORE ASSIGNMENT PHASE 

In order to develop a multiple linear regression equation 
for specifying the probability of fog occurrence, it was 
necessary to assign a quasi-continuous numerical value to 
each unique combination of ww, W, and CL in the synoptic re- 
port. The nature of the FOGCAT table facilitates this value 
assignment, as vell as the assignment of discrete, numerical 
or yes/no designations. 

The continuous value assignments to the FOGCAT elements 
are illustrated in Table V. The letter portion of each sub- 


category symool in Table IV was used to position the ww, W, 
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7 1:5 binatiıons into the listing order of Table V. For 
example, the four subcategories, bl thru b4, are combined into 
Mao This grouping of the four combinations that comprise 
subcategories bl thru b4 effectively gives each of these 

four combinations the same degree of certainty. 

The subcategory letter is paired with the visibility (VIS) 
value of each observation and a corresponding Base Score is 
determined. For example, an observation that has a VIS value 
of 92, and which receives a subcategory placement of kl or k2, 
would be assigned a Base Score of 360. As such, Table V was 
used as a decision logic table. 

It can be seen that a particular Base Score is assigned 
to many unique combinations of observed values of ww, W, CL, 
and VIS. For example, the Base Score of 260 can be assigned 
to any of 15 such combinations (five subcategories times the 
three visibility groups). The non-usage of the Base Scores 
of 160 and 180 makes a distinct gap or buffer zone in the 
continuous value range between those scores assigned to non- 
foggers and the ones assigned to observations with any 
possibility of being a fogger. 

The assigned Base Score was then scaled by dividing it 
by 5.8. This yielded the Fog Score, which was used as the 
predictand (dependent variable) during the regression phase. 
The Fog Score values (rounded) are listed next to the Base 


More values in Table V. 
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D. FURTHER REMARKS ON FOGCAT 

Other researchers would undoubtedly use different place- 
ments of the combinations of ww, W, and CL in the listing 
order of Table IV. Still others might suggest that the CL 
values could be used only if the sky cover code (N) value is 
co-evaluated. The Base Score assignment procedure could be 
kept "pure" by having a Base Score for each of the 72 combi- 
nations listed in Table IV and restructured in Table V. The 
Base Score range could then be 0-770 (intervals of 10 rather 
than 20), and/or the intervals could be variable so that Base 
Scores for foggers are close together but those for non- 
foggers would be farther apart (given consecutively listed 
subcategories). But, as a means to accomplish what is 
basically a feasibility study, FOGCAT is felt to be a "reason- 
able" scheme. 

Support for Nelson's (1972) general comments on the 
accuracy of the transient ship synoptic observations is shown 
in the counts listed in Table IV. According to the observing 
rules (U. S. Department of Commerce, 1969), those observations 
having present weather values in the 40's should also have VIS 
values in the 90-93 range. Subcategory a2 ended up with a 
total case count of 466. By the rules, all 466 of these ob- 
servations should have ended up being listed in the "Poor VIS" 
count for the a2 combination. Only 360 were so listed. It 
7014 have been possible to "bogus" in the proper visibility 


for those 106 faulty observations, but this vas not attempted. 
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The number of observations in Table IV for the a2 sub- 
category and the N major category have already been discussed. 
At least one other such count deserves mention, namely, the 
count of 313 observations for subcategory nl. The meaning 
o£ this unexpected high count is unclear, but the CL value 
of 6 suggests these observations may have been taken at or 
near fog regimes. A shift in the position in the listing 
order of Table IV for this subcategory should definitely have 


a noticeable effect on the resultant regression equations. 
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V. PREDICTOR INTERPOLATION 


The interpolation scheme chosen to provide values of the 
38 predictors (21 MOP's, 16 derived MOP's, and CLIMO) to the 
observation points is a natural bicubic spline curvilinear 
interpolation scheme, locally called SPLIN. The locally de- 
veloped computer program is available at the NPS W. R. Church 
Computer Center. The scheme can handle large arrays (up to 
300-by-300). But for time considerations, it was decided to 
limit the array to a 4-by-4 matrix for each interpolation 
performed. 

In order to use any interpolation scheme, it is necessary 
to be able to find the observation points in the same frame 
of reference as the predictors. This was accomplished through 
the latitude/longitude-to-I,J-coordinates conversion formulae 
described in Appendix F. The resultant I,J coordinates 
(63-by-63 grid) for a given observation point were used to 
generate the coordinates required with respect to the 249-by- 
249 grid used to derive the CLIMO values. 

The nature of the interpolation scheme, which is essen- 
PE Ti a Cubic curve fitting process, is such that errors can 
be introduced when the majority of the 4-by-4 array values are 
very similar to one another in sign and magnitude. These 


errors, and their possible effects, will be discussed later. 
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VI. REGRESSION SCHEME 


A. BACKGROUND 

A stepwise multiple linear regression program, called 
BMDP2R (University of California, 1973), is available at the 
NPS W. R. Church Computer Center, and was chosen as the means 
to derive a marine fog diagnostic equation. BMDP2R computes 
a sequence of multiple linear equations in a stepwise manner. 
At each step one variable is added to or deleted from the 
previous step's equation as dependent on the F-to-enter and 
F-to-remove criteria. Generally, the variable entered into 
the equation, in any given step, is the one which makes the 
greatest addition to the amount of predictand variance ex- 
plained (AR”). 

The BMDP2R program permits the specification as to vheth- 
er or not the y-intercept (YI) is to be treated as a variable 
or predictor. It is also possible to specify the tolerance 
level, which is a factor controlling the possible linear in- 
terdependencies or correlations between the various predictor 
variables. This tolerance specification is used by the pro- 
gram to inhibit the entering of a variable which has a high 
correlation to another previously entered variable. A case 
in point may be seen in Table VI, which lists correlations of 
all predictors to Fog Score along with the ranges of the 
variables. It can be seen that the correlations between the 


Fog Score predictand and both the EHF and SEHF predictors are 
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similar in magnitude and sign. The first regression equation 
obtained (see Table VII) used EHF as a predictor, but not 
SEHF. The correlation coefficient betveen EHF and SEHF is 
0.94 for the data used, vhich shovs the high, and expected 
1770 nd cenbatvesn these tvo predictor variables, and 
thence, vhy the tolerance control vill tend to keep SEHF out 
of the equation. The other correlations betveen the predic- 
tor variables are not shovn. 

The approach taken in this study vas to allov the treat- 
ment of the y-intercept (YI) as a variable and to use the 
BMDP2R program's default values of 4.0, 3.9, and 0.01 as the 
values for the minimum acceptable F-to-enter, maximum accept- 
able F-to-remove, and the minimum acceptable tolerance 


criteria, respectively. 


527 CULPUT RESULTS 
m Correlations 

The simple correlation coefficients (cc), between Fog 
Score and the 38 predictor variables, for the entire set of 
4481 observations are listed in Table VI. The groupings of 
the variables are for easier comparison between variables 
with similar characteristics, such as temperature, and/or 
model source. See Appendices C and D for model sources. 

The largest absolute cc value is the one listed for 
the EHF predictor. All of the other heat-flux type variables, 
except SOLARAD, also have relatively high cc values. 

Also vorth noting are the ccs for the FTR and CLIMO 


predictor variables. This leads to the expectation 
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that the results in Section VII (Verification Scoring) vould 
show that FTER possesses little, if any, skill over CLIMO in 
the diagnosis of the July 1976 marine fog regimes. 

2. Regression Equations 

The FTER and DDWW variables were removed from the 
BMDP2R program's use list before the generation of tne first 
set of regression equations was developed. FTER was removed 
because it is the existent predictor over which improvement 
is desired. DDWW was removed because it is a cyclic (0 = 360 
degrees) variable. Cycle variables are ill-conditioned for 
use in regression programs. Other cyclic variables, which 
were not recognized as such until late in the study, are PDW 
and SDW. Their inclusion as variables in some of the resul- 
tant equations may be undesirable, but time restrictions pro- 
hibited regeneration of the equations without them. 

The set of equations generated by BMDP2R for the 
whole study area (4481 cases) is listed in Table VII. They 
are listed in their respective stepwise order along with the 
variable regression coefficients, YI values, explanations of 


2 at the end of each step. 


variance (RZ), and the increase inR 
A discussion of Eq. (9.9) will occur later in the verification 
section. 

It is readily seen that the EHF variable, in conjunc- 
tion with YI, is the most significant predictor of the Fog 
Score. Since SEHF and the rest of the heat flux parameters 


were excluded by the tolerance specification, a separate run 


was made with EHF removed from the program's use list, with 
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the purpose of determining vhether SEHF vould be entered 
first in place of EHF. This vas expected, because SEHF has 
the second largest cc magnitude. The results (not shown) 
upheld the expectations. SEHF did replace EHF as the first 
variable to enter. 

Along with the EHF variable, only CLIMO, EAIR, and 
TSEA are significant variables in the explanation of the Fog 
Score variance. This significance is based on an R” increase 
Of 0.005 rather than the videly held standard of 0.01. 1f 
0.01 vere used as the significance limiter, EHF vould be the 
only such variable, along with the constant YI. 

None of the variables listed for Eq. (9) vere unex- 
pected. All of them can be logically and physically associa- 
ted vith the occurrence of marine fog. The association of 
T925 and fog is rather difficult to understand. This corre- 
lation is probably primarily due to the meridional distribu- 
tions of fog regimes and air temperatures, as is knovn to be 
the case vith TSEA. 

Tvo other equations generated for the vhole area are 
listed in Table VIII. The final step equation and its respec- 
tive coefficients are listed, as vell as the AR” due to each 
variable entered, the final R”, variables deleted (if any) 


Cheir associated AR” 


at entry and removal. Variables not 
given consideration (not in use list) during the regression 
are also listed. 


The development of Eq. (9T) stems from a transforma- 


tion of the Fog Score variable (predictand). This 
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transformation (T) takes the form of Fog Score = Fog Score - 
CLIMO. Effectively, this forces the regression coefficient 
of CLIMO to be one (1.0), as is reflected in Table VIII. 
Thus, CLIMO was not allowed entry into the equation as a pre- 
dictor during the stepwise procedure. The final ES E OE q . 
(9T) indicates that this approach holds little promise for 
useful application. 

The second equation listed in Table VIII is Eq. (9NC). 
In this experiment CLIMO was removed as a predictor variable. 
This regression result was run at the request of FNWC due to 
unavailability of acceptable fog frequency climatology for 
all oceans. The results, except for entry order, are similar 
to Eq. (9). The PDW variable is higher on the selection list 
than its counterpart in Eq. (9) and its significance is 
higher. In any case, the influence of this cyclic variable 
1s minimal. | 

Two sub-area equations were also obtained. The whole 
area file was sorted and copied into two smaller files, one 
(for Eq. (39)] for those cases (1504 observations) south of 
40N, and the other [for Eq. (40)] for all the other cases 
(2977 observations), namely, at or poleward of 40N. The final- 
step forms of these two equations are listed in Table IX. The 
final R” for Eq. (40) vas the highest obtained, vhile that for 
Eq. (39) was the lowest, with the exception of Eq. (9T). This 
“poor” showing by Eq. (39) might be attributed to the low 
occurrence of fog in the southern zonal belt below 40N. There 
were only 107 foggers (major category S and F) in the 1504 


observation reports in that area. 
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Two other runs made with BMDP2R are also worth men- 
tioning. One run had FTER as the only predictor variable in 
the use list, the other run had CLIMO as such. The runs were 
made against the whole area data file. The final R” for the 
BBB run is 0.5092, and tor the CLIMO run, 0.5391. This 
strengthened the expectation that FTER would possess little, 


if any, skill over CLIMO as a fog predictor. 
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VII. VERIFICATION SCORING 


A. BACKGROUND 

Two types of verification scores are used in the study to 
eet the “skill” of different predictors in describing the 
distribution of the FOGCAT major category designators (S, F, 
P, M, N). The first type of scoring is used in judging the 
effectiveness in forecasting ("novcasting" in this study) 
discrete (event, non-event) occurrences of a given phenomenon. 
The second type of scoring is used to judge the effectiveness 
of probabilities in estimating the occurrence of the event. 
In all, three scores are used in this study. The first two, 
Dën Skill Score (HSS) and Threat Score (TS), are of the 
first type of scoring approach. The third score, Probability 
Score (PS), evaluates the probability of fog as assigned by 
Fog Score, FTER, or CLIMO, and is of the second type described 
above. 

The exact formulae used for the three scores (HSS, TS, and 
PS) are given in Appendix G. The Heidke Skill Score has been 
used for quite some time in the meteorological community but 
may be falling into disfavor currently, while the Threat Score 
is currently in vogue with the National Weather Service. The 
Probability Score is appropriate to a predictor given in terms 


of a probability with a range of 0-100 per cent. 
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57 RESULT5 
1. Heidke Skill and Threat Scores 

Tables X through XIII show the various results for the 
HSS and TS calculations. The FOGCAT major categories of S 
and F were used as the criteria for classifying an observa- 
ión as a "fogger"”. The cutoff or cut value for a predictor 
is considered to be that value at which the estimate of the 
event is changed from a "yes event" estimate to a "no event" 
estimate (or vice versa). Only those cut values listed for 
FTER and CLIMO were used, because of the results of an 
earlier study (Renard, 1975b). The cut values for the equa- 
tions are those giving the best verification results. 

Table X reflects the results hopefully expected at 
the beginning of the study, namely, Equation (9)'s relatively 
Significant improvement over either FTER or CLIMO as a marine 
fog diagnostic tool. CLIMO and FTER are quite comparable 
in skill, which upholds the earlier expectation. Eq. (9T), 
which shows "poor" R” results, is comparable to Eq. (9) in 
ASS and TS results, as are Eqs. (9NC) and (9.9). The listing 
of the scores for Eqs. (1) through (8) shows the generally 
step-by-step increase in the scores. This listing complies 
with a suggestion by Panofsky and Brier (1968) to test each 
of the equations in the stepwise regression scheme. 

Listed in Table XI are the number of observations per 
day as well as the associated skill scores for the three pre- 
22 “ors FIER, CLIMO, and Eq. (9), using the best cutoff 


values from Table X for each of them. The day-by-day 
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comparison shows a good improvement of Eq.(9) over the capa- 
bilities of FTER and CLIMO in the diagnosing of the event. 
The day-by-day comparisons in Table XII for Threat Scores 
further highlights the merits of Eq. (9) in comparison to 
555R and CLIMO. 

The results listed in Table XIII for the tvo zonal 
belts again shov relatively significant improvement by the 
regression equations over FTER and CLIMO. The equations 
are comparable vithin the respective areas. That is, Eq. (9) 
(vhole area equation) exhibits about as much skill as the 
equations generated specifically for the sub-areas. 

Probability Scores 

These scores for the three predictor types are listed 
in Tables XIV and XV. CLIMO shows a slight improvement over 
FTER (the smaller the value, the better the skill) while 
Eqs. (9), (9T), and (9NC) show very little improvement over 
CLIMO. Eq. (9.9) shows the best improvement. This equation 
was obtained by changing the YI value of Eq. (9) with the 
express goal of optimizing the resultant PS. However, the 
optimization effect did not carry over to the HSS and TS 
results listed in Table X, 

The zonal-belt PS results listed in Table XIV show 
that the equations have a slight improvement in the northern 
belt over FTER and CLIMO. But in the southern belt CLIMO has 
the best score. This is probably due to the preponderance of 


“good'"' weather (non-fogger) observations south of 40N. 
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The daily PS values listed in Table XV support the 
findings in Table XIV, namely, that Eq. (9) does not shov 
much improvement over CLIMO and/or FTER. The daily perform- 
ance of Eq. (9.9) (optimized) shovs a far better than average 
improvement over Eq. (9), FTER, and CLIMO on a temporal fre- 
quency basis. 

The relatively small improvement of Eq. (9) over CLIMO 
and FTER in the PS results, as compared to the TS and HSS 
results, can be examined more closely by using Table XVI. 

Eq. (9) does a better job of specifying the foggers in the 
higher probability classes than CLIMO or FTER while CLIMO 

and FTER perform better in indicating non-foggers in the 

lower probability classes. The results of the PS calculations 
indicate these two forecasting aspects yield comparable PS 
values for CLIMO, FTER, and the regression equations. The 
figures listed in Table XVI compare favorably to those of 
Renard (19/5b) for July 1974 data, 

Whenever a predictor probability value was negative, 
zero was substituted so as not to cause complications with 
the PS formula. (See Table XIV.) This adjustment was not 


necessary for the HSS and TS calculations. 


30 





VIII. COMMENTS 


In section V it was indicated that the interpolation 
scheme might introduce errors into the data base. This error 
production can be readily seen by comparing the value ranges 
of the predictor parameters before and after the interpola- 
tion. The "before" values are listed in Appendices C and D, 
while the "after" values are listed in Table VI. The effect 
of this error introduction is clearly shown by the cases of 
below-zero value for the FTER variable listed in Table XIV. 
The FTER range has been expanded from Ü.Ü to -0.09 on the lov 
end and from 1.00 to 1.02 on the high end. The values may 
not seem significant, but 1183 observations vere assigned 
negative values of FTER by the interpolation process rather 
than the more desirable zero or small positive values. To 
check the extent of this effect, two additional Probability 
Score runs were made in which each of the 4481 values of 
FTER was modified. The results of these runs are the FTER 
+0.02 and FTER +0.05 entries in Table XIV. It can be seen 
that 891 of the 1183 observations received interpolated FTER 
values of 0.0 to -0.02, which is judged to be an acceptable 
deviation from zero. 

Mappings of sea-level pressure, EHF, FTER, probabilities 
from Eq. (9), and verifying observations by FOGCAT are dis- 
played for days 12 and 24 in Figures 5 through 14. The EHF 


charts are included since this model output parameter is 
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predominant in explaining a large amount of the variance of 
fog. Day 12 vas chosen because the scores produced for FTER 
and Eq. (9) are about "average", considering all days 
studied, vhile the performances of these tvo predictors for 


day 24 are "well above average". 
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IX. CONCLUSIONS AND RECOMMENDATIONS 


The primary objective of this study was to utilize the 
Model Output Statistics approach in the development of a sta- 
tistical regression equation that can accurately specify the 
distribution of marine fog. Through the interpretation, 
categorization, and statistical processing of various data 
elements, an equation was obtained whose variables are 
logically and/or physically related to fog (Table VII). The 
secondary objective, that of evaluating the skill of the 
different predictors (FTER, CLIMO, and the equations obtained) 
against a given data base, was fully achieved. 

The methodology of using diagnostic model output para- 
meters, rather than observational values, for describing the 
initial distribution of marine fog has tremendous potential. 
The equations obtained show significant increases in skill in 
performing this task in comparison to the use of climatology 
or the current FNWC fog forecasting product, FTER. FTER was 
Shown to have no appreciable skill over the NPS fog frequency 
climatology (CLIMO) for the criteria used. 

Regression Eq. (9), developed by using the whole area 
data file, proved to be as skillful in zonal band application 
as were the equations generated for the zonal bands. This 
Suggests that the zonal band partitioning does not segregate 


dissimilar areas well enough with respect to fog occurrence. 
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The following recommendations are offered for both short 


range and long range continuations of this work: 


(1) 


(2) 


(3) 


(4) 


5) 


(6) 


(7) 


(8) 


(9) 


Process the July 1976 parameters and observations 
again, using one or more different interpolation 
schemes. 

Break up the area of study into a multiple number of 
sub-areas on both geographical/meteorological bases 
and derive regression equations for each area. 
Recover the long wave (LW) portion of the THF para- 
meter and redo the July 1976 investigation including 
this parameter. 

Test the regression scheme in a forecast mode, using 
primitive equation prognostic model output parameters. 
Perform the same procedures, and modifications listed 
in these recommendations, on data from different 
months and/or synoptic times. 

Process the data bases using a discriminate analysis 
procedure. 

Analyze scatter diagrams from the data base using a 
non-linear approach. 

Plot the FOGCAT designations for each synoptic 

report each day and research the areas of fog 
generation, maintenance, and advection. Relate 

these fields to changes in the patterns of signifi- 
cant variables such as EHF. 

Experiment with changes to the FOGCAT structure and 


Fog Score assignment method or intervals. 


40 





(10) 


a) 


(12) 


Use Eq. (9) or (9NC) on a test operational basis as 
soon as practicable, and test the skill against FTER. 
Introduce a fog persistence parameter into any re- 
gression scheme for prognostic probabilities. This 
parameter also may be useful for the diagnostic 
probability equations. 

Eliminate the advection terms from consideration 
during the regression runs and redo the July 1976 
investigation. This would decrease the complexity 


of the parameters required for the "final" product. 
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APPENDIX A - Description of FNWC's Fog Probability 
Forecast Program (U. S. Naval Weather 
Service Command, 1975). 


1. Program Development 


The forecast for Advection Fog probability is based 
on a system of multi-parameter tests, vhere each parameter 
is veighted on the basis of test results. The veights of 
all parameters are summed, divided by 5, and normalized to 
give a value of zero to one, representing probability. 0.l 
would be 10% probability; 0.9 would be 90%. 

Advection Fog has five tests: 

(1) When surface air temperatures are near 0 
there is a high probability of fog. (Based on S. Petterssen) 
Weights favoring fog are assigned in the + 5°C range from 
Dec. 

(2) When condensation is taking place, there is 
usually an effective fog condition. This test is done by 
using evaporation calculations from moisture parameters. If 
evaporation is negative or zero, a certain weight is assigned 
to allov for the condensation factor. 

(3) If the near surface relative humidity exceeds a 
certain critical value (> 95%), a high weight is assigned on 
mae probability. 

(4) An air mass discriminant is applied. If the 
Surface geostrophic wind is being advected toward warmer 


(or colder) surface conditions, the air mass would be 
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classified accordingly--as in the old k and w air-mass system. 
A warm air mass would call for a high weight. The weight 
assigned is based upon calculations of the air mass wind com- 
ponent normal to the isotherms. 

(5) A veight is derived from application of a Fog 
Criterion. The criterion is based upon the change in dev 
point that occurs during the mixing of tvo air masses of dis- 
similar character. The resulting dev point depression in the 
mixed air is related to fog probability, vhich establishes a 
certain criterion value, a strong probability weight is 
assigned. 

User activities receive the fog probability forecasts 
in the form of a field, where each grid point has a fog 
probability value. 

2. Program Input 

(1) Sea-surface temperature analysis. 

(2) Surface air temperature analysis. 

(3) Surface vapor pressure analysis. 

(4) Sea-level pressure 6-hourly analysis. 

(5) P.E. prognoses of the above parameters from 0-48 
hours. 

3. Program Computation 

The various tests are calculated as described above 
for parameters at Tau's equal to 0, 6, 12 to 72 hours. The 
fog probability FTER is calculated for all land and sea 


points. 


43 





4, Program Limitations 


The five-parameter test method is admittedly some- 
what gross in character, but has enough built-in controls 
to insure a correct answer--given correct humidity values. 
The program gives excellent results in cases of broad-scale 
fog. The weakness of the program stems almost entirely 
from use of the large 63x63 grid. Local effects just cannot 
be accounted for, and fog is particularly susceptible to 
small-scale effects such as cold air drainage, small bodies 
of water, and man-made air pollution. Forecasters, in using 
the Advection Fog products, should alvays make allovances for 


50o0vn Local effects. 
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APPENDIX B: 


Station Number 


Ze . o ° 
eati on not Listed in reference. 


* 


32174 
5166 
717 
57707 
5215 


* 
82215 


32217 
82539 
32594 
32618 
70454 


kun 0 Used in study (U. S. Air Force, 


1772): 


45. 
46. 
46. 
48. 
50. 
20: 
Im 
ə. 
sır 
54. 
5 


oO UF ON OW ONY NN 
———- 


= Za 


Latitude 


“x 


Longitude Surface Elevation 
Cit) 

147.9 E 38 
Пун E 13 
Hp 26 
155 E 55 
577: 42 
50 2 TE unknown 
155.4 E Pr 
160.0 E 88 
SO sE 6 
166.2 E 19 
176. 7 W 14 


are as given with each synoptic 


Lat/long values listed 
report furnished by FNWC. 


“Location given slightly removed from actual location of 


55.2 N latitude, 166.0 E longitude to accommodate study 


area. 
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APPENDIX C. Output Parameter Descriptions. 


Source Model 


A. 


B. 


Symbol 


: = x Je 
Name Mean Minimm Maximum Units 


Description 
Analysis Parameters (FNWC's Mass Structure Model) 


PS 


TATR 


T925 


Sea Level Pressure: 1017.2 981.2 1034.0 (mb) 
Analysis of observed sea level parameter. 

Surface Air 

Temperature nx. 20 30.3 (56) 


Analysis of observation-level air temperature. 


Surface Vapor 

Pressure: jo 1 20 38.4 (mb) 
Analysis of observation-level vapor pressure derived 
from the dev point. 


925 mb Air 
Temperature: 1o72 -4.4 26.4 (2G) 
Analysis of 925 mb air temperature. 


Sea Surface 
Temperature: 16.6 2 29.0 Gc) 
Once-daily analysis of observed sea-surface temperature. 


P.E. Parameters (FNWC's Primitive Equation Model) (Kaitala, 1974) 


dE 


SOLARAD 


Surface Air 

Temperature: 19,5 . 28.2 Ce) 
Derived from surface air and potential temperatures, 
boundary layer depth, upper-level winds extrapolated to 
surface, air density, drag coefficient, gustiness factor, 
and empirical constants. 


Surface Vapor 
Pressure: Tee 6.6 927 (mb) 
Derived from model"s mixing ratio. 


Solar Radiation: 405 25.9 88.4 5: 

cm? +H) 
Calculated absorption of incoming short-wave (solar) 
radiation (positive downward) (H = depth unit). 


Evaporative Heat 

Flux: 3.5 -14.9 25 Lë: "CH 
cm! -H) 

Derived using air density, drag coefficient, extrapolated 

winds, and mixing ratios. 
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APPENDIX C (continued) 


SEHF 


Sensible/Evaporative 

Heat Flux: MONO 277 Ge) 
cm?» H) 

SEHF = Sensible Heat Flux (SHF) + EHF 

Sensible Heat Flux: -5.6 -28.4 3.6 Lë, "CH 
cm? +H) 


Recovered from SHF = SEHF-EHF. Originally derived by FNWC 
using drag coefficient, extrapolated winds, surface air 
temperature, TX, density, and constants. 


Total Heat Flux: -37.6 -80.9 -9.0 Lë: "CH 
cm” +H) 

THF = SEHF - SOLARAD + LW, where LW is the heating due 

to long-wave (terrestrial) radiation. 


C. Marine Wind Model (FNWC) 


VVWW 
DDWW 


Marine Wind Speed: 1953 0.0 74.8 (knots) 

Marine Wind 

Direction: 18.7 0.0 36.0 (degrees/ 
10) 


Both variables derived from a dynamic balancing of surface 
wind and sea-level pressure. 


D. Spectral Ocean Wave Model (5.0.W.M.) (FNWC) 


HW 


PPW 
PDW 


WCP 


Significant Wave 


Height: 4.4 0.0 2 (feet) 

Primary Wave Period: 2 0.0 20.0 (sec) 

Primary Wave 

Direction: 12 0.0 57 (degrees / 
10) 

secondary Wave 

Period: 5.6 0.0 25.0 (sec) 

Secondary Wave 

Direction 11,5 050 57 (degrees/ 
10) 

Probability of 

White Caps: 0.3 0.0 7.0 (per 
cent) 
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APPENDIX C (continued) 


E. Other Model Output Parameters (FNWC) 


SSTA 


Sea Surface Tempera- 

ture Anomaly: -0.6 -4.5 77 eo) 
Calculated anomaly of sea-surface temperature from the 
mean of the day as interpolated from the monthly mean 


values. 


Probability of Fog: 0:1 0.0 EO (per cent) 
Forecast of advection fog probability. 


“These values are for the grid points in the rectangular area (Fig. 3) 
encompassing the area of study for 1-30 July 1976, including 0.0 for 
wave parameters over land points. 
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Units 


(m/sec) 


(m/sec) 


(m/sec) 


(m/sec) 


(°K) 


(°K) 


(°K/mb) 


(°K/mb) 


(°C) 


(°C) 


5 
Hour) 


APPENDIX D: Derived Parameter Descriptions 
* * 
Symbol Name Mean Minimm Maximm 
Description 
u Zonal Wind 
Component 1.7 -20.9 26.9 
= -VVW sin (DOW - 10). 
V Meridional Wind 
Component fo eal 25.4 
v = -VVW cos (DDWW + 10). 
CAPU I Directional 
Vind Component: 12 7222485 2055 
CAPU = -u » sin(LNGA) - v = cos(LNGA), (Haltiner, 1971). 
CAPV J Directional 
Wind Component: -2.0 -35.1 Hö 
CAPV = u - cos(LNGA) - v “ sin(1NGA), (Haltiner, 1971). 
where LNGA = -10 -(I,J point longitude).”” 
THETAX Potential 
Temperature X: 2 000775 30077 
Derived using PS, TX. 
THETAR Potential 
Temperature R: 289.5 276.5 7 
Derİved using PS, TAIR. 
STABX Stability X: -0.00 -0.19 0.07 
Derived using [THETAX -(THETA of T925)]/(PS-925). 
Value greater than zero indicates absolute instability. 
STABR Stability R: -0.07 -0.25 Om 
Derived using [THETAR -(THETA of T925)]/(PS-925). 
Same value effect as STABX. 
ASTDX Air-Sea Temperature 
Difference X: 1.9 -3.4 AS 
ASTDX = TX - TSEA. 
ASTDR Air-Sea Temperature 
Difference R: kə -5.6 27 
ASTDR = TAIR - TSEA. 
ADTSEA Advection of TSEA: 0.04 -0.79 0 95 
Formulae and notes below. 
ADTX Advection of TX: 10 1.02 


Formilae and notes below. 
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APPENDIX D (continued) 
anes 


ADTATR Advection of TATR: 0.04 -0.85 0.86 (567 
Hour) 
Formulae and notes below. 
y 
AASTDX Advection of ASTDX: -0.00 -0.69 0.78 (Ser. 
Hour). 
Formilae and notes below. 
sock 
AASTDR Advection of ASTDR: 0.00 -1.02 0.64 (2G; 
Hour) 


Formulae and notes below. 


“These values are for the grid points in the rectangular area (see Fig. 3) 
encompassing the area of study. 


tek 
I,J point longitude calculated using I,J-to-latitude/longitude 
conversion formulae described in Appendix E. 


IN 


Mean may be biased towards zero. Only values used in study were 
calculated. Grid points outside study area received zero. 


Advection Formulae and Conditions: 


For the advection of a quantity (R) the formula, ADQ = -V-V(Q), was 
used in the folloving finite difference form: 


RMAP 


ADQ = – — Бӱ (САРЈ" (Оз - 91.1) T+ GAP: Ui +87 129)» 
27277:5:5 0 07 ein Clatitude)) 
səl 555 0 a 7 “ın (60) (31205 


(31.205 = grid mesh lengths, pole to equator, on FNWC's I,J grid). 


In the temperature advection calculation for point C, using the five 

grid points illustrated below, one or two of the points, namely, A, B, 

D, or E, may be outside the study area. In the bogusing method 

suggested by Mr. Leo Clarke, FNWC, when a non-center point (e.g., point A) 
was judged to probably produce a land and/or dissimilar sea area 
influence on the resulting advection, the center point value (point C) 
was substituted for it. This 'bogusing'' ls necessary to maintain a 
"purely" marine characteristic in the resultant parameter value. 


The set of study area boundary grid points used for bogusing (and 
some double bogusing) are depicted in Fig. 2. The upper case letters 
denote the points whose values were used for the bogusing. The lower 
case letters mark the positions of the adjacent points being bogused. 
The only boundary point close enough to land to give concern is to 
the one near station number 32594 (south tip, Kamchatka Peninsula). 


50 





APPENDIX D (continued) 


However, it was judged to be over the ocean after examining TX, TAIR, 
and TSEA values at this and adjacent land/ocean I,J-grid points. 


A 
HER 
B С D 
GE 1,3 Dol 
I 
Tə 
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APPENDIX E: Conversion of I,J grid coordinates to latitude/ 
longitude positional values (FNWC Subroutine 
Library). 


The following formulae were programmed in Fortran as a 
subroutine used by several different large programs developed 
mer this study: 


LONGITUDE = K - arctan 167-97-27 (1-13) 


5:2 07485 (JJ 2 
LATITUDE = arcsin 7 7 7 
(RED) “° + (1-1 + 5: 


where: K = -10 or 1/0 depending on quadrant used in 
study. (-longitude means West longitude) 
J_ = Pole position = au 
* 
I. = Pole position = 31 


RED = Distance from pole to equator in grid mesh 
lengths = 31.205 


* 
32 used on IBM 360/70 at NPS W. R. Church 
Computer Center. 
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DSEPENDIX E: 


Conversion of latitude/longitude position to I,J 


grid coordinates (FNWC Subroutine Library). 


The following expressions were programmed for use during 
the interpolation phase of the study: 


A : cos 6 : S 
I = L. + [RED та. Gos 510571 
= a Cos a0 | ә 
where $ = Latitude 
A = Longitude (negative = West) 
* 
I, = Pole position = 31 
J, = Pole position = 3 
RED = Distance from pole to equator in grid mesh 


lengths = 


1502 


7 used on IBM 360/70 at HPS W. R. Church 


Computer Center. 
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APPENDIX G: Verification Score Formulae 


From a contingency table and associated information the 
Heidke Skill Score (HSS) and Threat Score (TS) can be 
calculated: 


EVENT ESTIMATED 


ES NO 

YES A C Total (T) = A + B + C + D 
EVENT 
OBSERVED Ho OE Correct 

NO B D Forecast (FC) = A + D 

wee EX 52585 m: 
HSS = runs x Range: m xn < 56 > 1 
B+ C 

where EX = R RE ee : No. of expected 


correct foracasts due to chance. 


A A 


ТӘ = 725” BC 


Both scores indicate more skill with larger 
positive values. 


The Probability Score (PS) is from that given by 
Panofsky and Brier (1958) and may be written as (Renard, 1975): 


no ny 
PS =- HN 9 +) 0-29 
i=1 j=1 
Ranges SS SZ 
where N = Total mumber of cases 
n. = Total number of non-events 


p; = Associated probability value for the non-event 
пј = Total number of events 
p. = Associated probability value for the event 


The closer to zero, the greater the skill. 
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TABLE I. 


Abridged version of internationally used veather-code figures 
and definitions for reporting present and past veather and 
lov clouds in the surface synoptic report (U. S. Departments 
of Commerce, Defense, and Transportation, 1969). 


Present Weather 


Code 
Value 


00-03 


04-09 
10 


11-12 
13-17 


18-19 


20 
21-23 


24 
25-27 
28 
29 


Definition 


Characteristic change 
of the state of the sky 
(cloud) during the past 
hour. 

Haze, dust, sand, or 
smoke. 

Deep light fog. 

Shallow heavy fog. 
Lightning, thunder, or 
precipitation within 
signt, not reaching the 
ground. 

Squall(s), fumel 
cloud(s) during the 
past hour. 

Drizzle during the 
past hour, 

Rain, snow, or rain and 
snow during the past 
hour. 

Freezing drizzle during 
the past hour, 
Shower(s) during the 
preceding hour. 

Fog during the past 
hour. 

Thunderstorm during the 
past hour. 


Present Weather 


Code 
Value 


30-39 
40 


41-49 


50-59 


60-63 
64-65 


67 
68 
69 
70-79 
90-99 


55 


Definition 


Duststorm, sandstorm, 
drifting or blowing snow. 
Fog at distance, but not 
at station during the 
past hour (visibility 
less than 1 km). 

Deep heavy fog at the 
time of observation 
(visibility less than 

1 km). 

Drizzle, or drizzle and 
rain. 

Slight to moderate rain. 
Heavy rain. 

Slight freezing rain. 
Moderate or heavy 
freezing rain, 

slight rain or drizzle 
and snow. 

Moderate or heavy rain 
or drizzle and snow. 
solid precipitation not 
in showers. 

Showery precipitation 
or precipitation with 
current or recent 
thunderstorms. 





TABIE I. (continued) 


Past Weather 


Code 

Value Definition 

0 Cloud covering % or less 
of sky throughout the 
period. 

L Cloud covering more than 
# of sky during part of 
period. 

2 Cloud covering more than 
% of sky throughout the 
period. 

3 Sandstorm, or duststom, 
or bloving snov. 

A Heavy fog, thick haze, 
or smoke. 

о Drizzle 

6 Rain 

7 Snow, rain and snow 
mixed, or ice pellets. 

8 Shower (s) 

9 Thumderstorm, with or 
without precipitation. 

Visibility 

90 Less than 50m 

91 0-199 m 

7 200-499 m 

> 500 m - 0.99 km 

94 [= 1239) lan 

25 2 - 3.99 km 

96 4 - 9,99 km 

9/-99 equal to or greater than 


LO km. 


Low Cloud Type 


Code 
Value 


0 


Q ND Un FP0-(o N 


NO 


56 


Definition 


No low clouds 

Ragged cumulus o£ fair 
weather. 

Generally towering 
cumulus , 

Cumulonimbus vithout 
cirriform or anvil tops. 
Stratocumulus formed by 
cumulus spreading out. 
Stratocumulus not formed 
by cumulus spreading. 
Stratus or fractostratus. 
Fractostratus of bad 
weather. 

Cumulus and stratocumulus, 
with bases at different 
levels. 

Cumulonimbus with 

emr Tusi GOD. 

Low cloud obscured. 





" (2)67-07 TO ‘gZ ‘ZI ‘TI ‘OT 0U sT (Mmm üəüM BOF AFTIUSPT 07 posn sepod AITTIGISTAry 
"“әлпатј әроә YLM 1584 до 4uəsəqd ÁQ P9IPOTPpur UOTIDTA1IS9A ÁITTIQISTA JO AÇI IV x 
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TABLE III. Groupings and symbols used in FOGCAT 
categorization scheme. See Table IV. 


Present Weather (ww) Past Weather (W) Low Cloud (CL) 

Symbol Associated Symbol Associated Symbol Associated 
ww Codes W Codes CL Codes 

41G 41-49 4 4 6 6 

LOG 10,28,40 4,5 4,5 27 Dw 

11G 11,12,20,24 5G O LS B / 

50G 50-59 2G 0 x any CU mot 


listed above 
60G 60-63,66,68 Bi any W not 
listed above 
* any ww not 
listed above 


G = Group. B = Low clouds obscurred. 
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TABLE IV. Scheme for categorizing observations according to 
likelihood of fog (FOGCAT). See Table III for 
symbols. Number of July 1976 Worth Pacific Ocean 
observations in each subcategory stratified by 
visibility code. 


Low Number of Visibility (VIS) 
Major Sub- Present Past Cloud Observations by VIS Code 
Category Category Weather Weather Type ATI Fair Poor 


(ww) (W) (CL) (90-99) (94-95) (90-93) 
Strong al &1G 4,2 6 66 19 by) 
Foggers a2 5 " B 466 83 360 
= S bl di 2G 6 24 7 2 
b2 o ` B 28 6 1A 
b3 10G 4,5 6 “1 9 2 
b4 iz 4 B 36 12 16 
el 41G x 6 2 L 1 
ez . 4,5 57 26 5 13 
c3 2 Ҹ̧ * 64 14 39 
c4 10G 2G 6 30 8 2 
Total S cases: — 
Foggers dl 41G * B 1 3 17 
= F d2 y 2G 27 12 2 3 
d3 50G 4 6 ИЛ 3 0 
dA S ú B 19 5 6 
də 10G Bi 6 3 2 0 
d6 a 2G B 13 4 1 
d7 É & Del 14 2 1 
е1 416 x R 5 2 2 
e2 : 2G * 16 4 5 
e3 50G 4 57 3 0 0 
ед 10G * B 0 1 
e5 2 & * 37 9 1 
e6 11G & 6 6 0 0 
e7 ğ s B 7 2 1 
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TABLE IV, 


zü 
f2 
f3 
£4 
Lə 
F gl 
g2 
g3 


ва 
gö 
hl 
h2 
r 
iz 


ıl 
Total F cases: 


Past/Weak ` kl 
Foggers k2 
= P “ 

22 


ml 
m2 
nl 
ol 
о2 
pl 
pe 


p3 
Total P cases: 


(continued) 


&1G 


60 


57 


m Ln Oh D C, ON Ww QON ON 


15 
12 


24 


D 


N N OM OF NY О (90 О 34 О 49 О О 


O O rəm © OC RESTA 


Go O O O O mə Ñ F: Ñ m o o 





TABLE IV. (continued) 


Maybe 
Foggers 
El 


ql 
q2 
q3 
q4 
q) 
q6 
q/ 
El 
иә 
E3 
rA 
ə 
r6 


Total H cases: 


Non-Foggers ul 


= H 


u2 
vl 
vl 
xl 
vl 
y2 
zı 


Total N cases: 


CO 


ə 
29 


ə Un 
Y N 
— 


ES 


14 
“ 


5 
19 
19 


Uu wa ooo OA O + 
y 


N 


244 


8 
E 
SS 
B 13 
Diy) Ад 
3 7 687 
5,7 1541 


DI x x 


* 102 


2604. 


_50 


H OU Cə 


rem вә KO (иә O 


Note: A dash in a visibility scheme indicates that particular count 


vas not done. 
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TABLE V. Major fog categories and sub-category fog groups 
with associated Base and Fog Scores as a function 
SiO: 


Visibility Code Values 


96-99 94-95 90-93 Base Fog 
(Good) (Fair) (Poor) Score Score 
Major 
Category Subcategory Groups 
a 580 100.0 
a b 560 50 
a b 2 540 5 
S b “se. d 2 89.7 
5 Zen d e 500 207 
d e in 480 625: 
e f g 460 EDS 
£ g h 440 5 
F g h 1 &20 72.4 
h il sə. 400 890 
d bt 360 6560 
5 65 k 360 6251 
k k i) 340 55 
2 g m SAY, 5057 
m m n 300 5177 
P n n O 280 | 48.3 
O O ити 260 44.8 
p Son 240 41.4 
= a q q 220 5779 
M E sr Ели 200 34.5 
кос 5 5 5 180 0 
used A ee ır 160 2770 
u 140 24.1 
u u V T20 2011 
V V W 100 7. 
VW W X 80 10 
N x x y 60 de 
y Z 40 6.9 
y 2. 20 3,4 
Z 0 20 
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TABLE VIII. Regression equation predictor parameters and coefficients, 
and associated explained predictand variances (R”) for 
special treatment of climatology parameter (final-step 
form), vhole area, North Pacific Ocean (about 30-60 N), 


July 1976. 
Equation 9T Equation 9NC 
Fog Score transformed (T) Regression Variable AR” 
(Fog Score = Fog Score - CLIMD) Coefficient (in entry 
before regression. order) 
Regression Variable AR” 5: A —.. 
Coefficient (in entry -1.8859 EHF 0.1183 
order) -4.4919 ASTDR 0.0035 

1. CLIMO - 2120620 ADTATR 0.0042 
-1.8328 EHF" 0.0299 -4.9522 TSEA 0.0028 
+2.9127 FAIR 0.0279 +2.9131 FAIR" 0.0091 
-2.0207 TAIR“ 0.0154 40.9988 T925 0.0028 
2075283 V O 0046  -0:2555 PIW 0.0012 
-1.3244 ASTDR 0.0022 +0.2698 V 0.0004 
- YI upon entry 0.0557 

- YI upon removal -0.0004 
Final R” - 0.1353 Final R” - 0.5843 
HER, DDWW, (SEHF, THF, SER, CLO, FIER, DDWW, (SEHF, THF, 
CAPU, CAPV) ^ not in the use SHF, CAPU, CAPV) ^“ not in the 
list. use list. 

Z 


* Most Significant Variables (based on AR” > 0.005). 
“* Omitted through error. 
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TABLE IX. 


Regression equation predictor parameters and coefficients, 
and associated explained predictand variances (R2), final- 
step form, for two sub-areas, North Pacific Ocem, July 1976. 


Equation 39 Equation 40 
Based on 1504 observations south Based on 2977 observations at and 
of 40 N. : north of 40 N. 
Regression Variable AR” Regression Variable AR” 
Coefficient (in entry Coefficient (in entry 
order) order) 
ҹа * 
-0.9826 EHF 0.0929 21.8849 ar 055577 
: * 
0.6595 CLIMO 0.0246 -2.0569 EHF 0.0601 
* “ls 
17 5560 EATR 0.0089 -7.3865 ASTDR 0.0058 
=! 355 SSTA 0.0025 0.2945 CLIMO 0.0020 
* 
-0.1608 SDW 0.0016 -7.6194 TSEA 0.0165 
-74.0931 STABR 0.0018 0.7898 T925 0.0014 
- YI upon entry 0.2842 0.5861 V 0.0011 
- YI upon removal -0.0003 -26.5505 AASTDX 0.0010 
HW upon entry 0.0043 ADTSEA | 
HW upon removal -0.0014 57 0.0027 
ADTSEA 
upon removal -0.0004 
Final R* = 0.4219 Bıı 0 00/68 
FIER, DDW, (SEHF, THF, SHF, 
CAPU, CAPV)”” not in use list 
for both equations. 
2 


* Most Significant Variables (based on AR” > 0.005). 


** Omitted 


through error. 


@ Not best equation (39) possible. A full use list mm (not shown) 
made late in the study gave SEHF a higher initial correlation to 
Fog Score. 
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TABLE X. Heidke Skill and Threat Scores, whole area, HNorth 
Pacific Ocean (about 30-60 N), July 1976. FOGCAT 
designators S and F used to identify marine fog 
observations. See Table IV and text. 


Cle Ost 
Predictor (Cut) Skill Score Threat Score 
(Per cent) 

FTER 30 iL 0.304 

70 029 -023 
CLIMO 30 771 -325 

40 6232 57271 
Egel 40 .340 2070 
EG pe 40 55 "572 
Bas 40 .346 0365 
Еа 4 40 550 77: 
Bass 40 701 77 
Eq 6 40 -3/6 -338 
Egay 40 .378 . 389 
Eq 8” 40 ‚386 ‚394 
Eq 9 40 75: 502 

45 55 — 
Eq 9.9 25 -3380 72 
Eq 9T 45 “2 2277 
Eq 9NC 40 79 . 384 


Scoring performed on vhole area data file, i.e., 4481 cases, 
of which 1055 are classified Foggers, using major categories 
S and F as criteria. 


+ 
25: . one varıable's coefficient makes these results 
slightly suspect. 
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TABLE XI. Daily Heidke Skill Scores, whole area, North 
Poesie Ocean (about 30 - 60 N), July 1976. 
FOGCAT designators S and F used to identify 
marine fog observations. See Table IV and text. 


DAY Cases Foggers ETER CLIMO Ед 1, 
Cut=30 =30 =40 
(Per cent) 
1 125 28 US 9.356 0.459 
2 162 34 ui .346 273 
3 145 85 .140 ‚245 82 
4 147 33 „182 121: .379 
5 152 39 OG „149 5775: 
6 156 30 -0.060 . 200 :7- 
7 148 7 . 027 ‚160 әр 
8 140 35 . 050 322 7 
> 160 32 -219 :326 . 408 
10 163 38 . 459 302 .540 
il 192 25 . 284 2209 319 
12 lr. 33 H ao .442 
13 T99 32 Zone 290 “215 
14 196 ди -OTO 17” 201 
ES 160 34 “225 , 386 (455 
16 156 34 „228 „328 -363 
Ly i: 31 2 . 200 336 
18 164 38 277 . 148 . 402 
19 169 40 . 414 .264 5550 
20 147 32 (450 - 308 ‚446 
21 138 24 “io 7 7: 
22 E 40 ‚407 ‚416 .438 
23 148 40 267 . 375 ‚483 
24 157 32 . 488 3806 2009 
25 143 45 ‚444 2507 1012 
26 156 42 „455 ‚490 ‚451 
27 146 41 75: 751 577 
28 143 44 .348 107 1777 
29 107 49 . 3095 212 .402 
30 259 49 26m 2161 5:7: 
Total 4481 1055 


CLIMO has more skill on 17 of the 30 days compared to FTER. 
Eq. (9) has more skill on 24 of the 30 days compared to FIER, 
Eq.(9) has more skill on 28 of the 30 days compared to CLIMO. 
Eq.(9) has more skill on 23 of the 30 days compared to CLIMO 
and FIER. 
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TABLE XII 


CLIMO has 
Eq. (9)has 
Eq. (9)has 
Eq. (9)has 


Daily Threat Scores, whole area, North Pacific 


Ocean (about 30-60 N), July 1976. FOGCAT 
designators S and F used to identify marine fog 
See Table IV and text. 


observations. 


DAY 


| 
O xO Oo — OYU fo ER 


PRA pa A 
Pp Mn ә 


— = 
Co Y OU 


more 
more 
more 
more 


FTER CLIMO 
Cut=30 -30 
(Per cent) 
0.241 0.346 
120. 2.01 
257 5: 
5- 7: 
55 ‚276 
„134 2297 
5156 ‚224 
. 194 , 349 
7: . 329 
2621 757: 
76 . 305 
MODO . 349 
5720 5 
55 KA 
36 . 379 
1520- „3935 
233 , 269 
. 265 250 
‚354 22 
¿212 . 340 
5170 —. 
. 400 426 
‚429 . 409 
‚464 5200 
77. 57 
„457 ‚486 
‚407 рз 
‚388 . 284 
.514 .347 
7: ‚304 
skill on 19 of the 30 days 
skill on 29 of the 30 days 
skill on 27 of the 30 days 
skill on 26 of the 30 days 
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EOI) 
=40 


0.446 
700 
. 342 
5505 
.314 
‚247 
ze 
. 350 
‚374 
. 491 
57 
. 421 
„241 
2.51 
. 431 
572 
, 362 
. 400 
. 49) 
. 424 
. 300 
. 439 
(477 
. 309 
75 
. 4530 
‚414 
‚410 
„451 
(426 


compared to FTER. 
compared to FTER. 
compared to CLIMO. 
compared to CLIMO 

and FTER. 





TABLE XIII. Heidke Skill and Threat Scores, sub-areas, North 
Pacific Ocean (about 30-60 N), July 1976. FOGCAT 
designators S and F used to identify marine fog 
observations. See Table IV and text. 


Predictor Cut Skill Score Threat Score 
(Per cent) 


Area: Northern Belt” (2977 observations, 948 foggers) 


FTER 30 0.233 0.335 
70 026 026 
CLIMO 30 151 341 
40 Ee 259 
Eq. (9) 40 299 412 
45 DD 397 
Eq. (40) 40 306 418 
45 57 412 


Area: Southern Belt” (1504 observations; 107 foggers) 


FTER 30 0.156 | 0.133 
70 -0.001 . 000 
CLIMO 30 .067 .043 
40 . 000 .000 
Еа. (9) 35 .288 .213 
40 ‚247 .174 
Eq. (39) 30 283 . 209 
35 . 254 илә 


* Northern Belt 
**Southern Belt = whole area cases south of 40 N 


whole area cases at or north of 40 N 
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TABLE XIV. Probability Scores (PS) for whole and sub-areas 
of the North Pacific Ocean (30-60 N), July 1976. 
FOGCAT designators S and F used to identify 
marine fog observations. See Table IV and text. 


Whole_Area 
No. of observations 

Eredictor PS Max Min vitü value below zero 
FTER 0.350 OZ -0.09 1183 

FTER +0.02 . 349 1.04 -0.07 292 

FTER +0.05 1027 1:07 -0.04 21 

CLIMO SES 20.0 0.0 0 

Bar (1) 5325 96. -A8. 195 

En (2) ЈӘИӘ 57: -36, T95 

Ea. (3) ¿317 87. -33, 185 

Eq. (4) 2.7 88. -35, 209 

Bar (5) wore 94. -37, 192 

Bar (6) 2515 35, -53. 257 

Bar (7) ome? 88 -53. 2:2 

Edu (8)* 127 88, -52. 268 

Bar (9) Sal 89. -50, 253 

Bar: (9.9) 290 77 -60, 576 

cee (9T) 0 50 -55. 261 

LG (9NC) la 77 -A47, 277 


Northern Belt (cases at or north of 40 N) 


Predictor PS Predictor PS 

FTER 0.446 el, 02.393 
кҝ CLİMO......0.412.......8q. (40) 0:3802 
Southern Belt (cases south of 40 N) 

FTER 07162 Pa m (9) 0.149 

CLUMO 0.124 EGRESO) OFS) 


* 
Error in coefficient makes result slightly suspect. 
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iH 57077 007 Po abiliEy Scores, whole area, North Pacific 
Ocean (about 30-60 N), July 1976. FOGCAT designa- 
tors S and F used to identify marine fog observa- 
tions. See Table IV and text, 


DAY FTER CLIMO Eq. (9) Eq. (9.9) 
1 0.384 o 231 0.291 oğ 
2 KEN „318 5. . 289 
3 . 394 „341 — ‚326 
4 „434 505 . 328 . 293 
5 „471 , 369 385 364 
6 447 , 299 ‚342 ‚246 
/ 77 5279 516 577- 
8 WIER .314 „332 ‚308 
9 EEN .276 (285 . 260 

10 „281 . 301 .284 5257 
ni .299 јә . 303 ‚262 
12 340 25091 , 298 „270 
13 „302 . 309 . 300 . 280 
14 ‚305 KE 750 „245 
15 50 . 265 211 .253 
16 1352 , 299 ‚316 288 
17 366 . 308 .320 . 296 
18 365 .328 , 298 , 288 
19 5) — ‚276 ‚265 
20 5525 ии. .305 .278 7 
21 5 . 282 .299 Dr, 
22 ‚334 2007 ASS . 297 
23 .337 .326 .292 .289 
24 .267 .370 75 577 
25 ли „351 , 308 .307 
26 .307 .297 .291 282 
27 ‚332 .366 .368 „341 
28 „369 ‚390 362 75: 
29 . 344 ‚394 “2339 5 
30 „450 .391 .348 .354 


CLIMO had more skill on 22 of the 30 days compared to FTER. 
Eq. (9) had more skill on 27 of the 30 days compared to FTER. 
Eq. (9) had more skill on 14 of the 30 days compared to CLIMO. 
Eq. (9) had more skill on 13 of the 30 days compared to CLIMO 
. and FTER. 
Eq.(9.9)had more skill on 29 of the 30 days compared to FTER. 
Eq.(9.9)had more skill on 30 of the 30 days compared to CLIMO. 
Eq. (9.9)had more skill on 27 of the 30 days compared to Eq. (9). 
Eq.(9.9)had more skill on 26 of the 30 days compared to FTER, 
CEMO Tan Eqo): 
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Figure 3. Fleet Numerical Weather Central's 63x63 grid, 
with outline of North Pacific Ocean rectangu- 
lar grid area used in study. See text. 
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